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It has been proposed that localized zero-energy Majorana states can be realized in a two- 
dimensional network of quasi-one-dimensional semiconductor wires that are proximity- coupled to a 
bulk superconductor. The wires should have strong spin-orbit coupling with appropriate symmetry, 
and their electrons should be partially polarized by a strong Zeeman field. Then, if the Fermi level 
is in an appropriate range, the wire can be in a topological superconducting phase, with Majorana 
states that occur at wire ends and at Y junctions, where three topological superconductor segments 
may be joined. Here we generalize these ideas to consider a three-dimensional network. The po- 
sitions of Majorana states can be manipulated, and their non-Abelian properties made visible, by 
using external gates to selectively deplete portions of the network, or by physically connecting and 
redividing wire segments. Majorana states can also be manipulated by reorientations of the Zeeman 
field on a wire segment, by physically rotating the wire about almost any axis, or by evolution of 
the phase of the order parameter in the proximity- coupled superconductor. We show how to keep 
track of sign changes in the zero-energy Hilbert space during adiabatic manipulations by monitoring 
the evolution of each Majorana state separately, rather than keeping track of the braiding of all 
possible pairs. This has conceptual advantages in the case of a three-dimensional network, and may 
be computationally useful even in two dimensions, if large numbers of Majorana sites are involved. 



I. INTRODUCTION 

The last several years have witnessed extensive re- 
search into the possibility of realizing zero-energy 
Majorana-fermion states ("Major anas") in condensed 
matter systems. Driven by the interest in non-abelian 
states of matter and by possible implementations for 
topological quantum computation^'^, several physical 
systems have been arenas for the search for Majorana 
fermions. Initially, these systems were primarily two 
dimensional, such as the u = 5/2 fractional quantum 
Hall state, in which Majorana fermions are believed 
to be attached to quarter-charged quasi-particles and 
quasi- holes in the PfafRan or anti-Pfaffian states^'^; two- 
dimensional spin polarized p-wave superconductors^'^; 
surface states of three-dimensional topological insula- 
tors coupled to superconductors^; vortex cores in hybrid 
systems made of semiconductors with strong spin-orbit 
interaction in proximity-coupling to superconductors^'^; 
and more. While most of the research has been theoreti- 
cal, interferometry experiments (according to suggestions 
such as Refs. 10 and 11) have yielded data that may be 
interpreted as resulting from the presence of such excita- 



tions-^^ 

One dimensional systems have been predicted to host 
Majorana modes as well. Ideas of end-point Majoranas 
in spin polarized p-wave superconducting wires were al- 
ready discussed by Kitaev^^ nearly a decade ago using 
a toy tight binding model for a spin-polarized p-wave 
superconductor. Possible realizations in semiconductor 
quantum wires with strong spin-orbit coupling, such as 
In As nanowires, in the presence of a strong Zeeman field 
and proximity-coupling to a bulk s-wave superconductor, 
and their manipulation by means of external gates po- 
tential, were suggested in Refs. 16-19. One-dimensional 
systems of the required form might also be realized in 
edge states formed at the surface of a topological insula- 
tor coupled to a superconductor and a ferromagnet^^'^^. 

The presence of localized zero energy Majorana 
fermions states separated from each other by distances 
large compared to the coherence length makes the ground 
state of the system degenerate, in all dimensions. It 
is well known that in two-dimensional systems braid- 
ing Majorana modes around each other produces unitary 
transformations within the ground state subspace which 
generally do not commute; hence the statistics is non- 



2 



Abelian^^'^^. More recently, it has been shown that sim- 
ilar manipulations leading to non-Abelian statistics are 
possible in two-dimensional networks synthesized from 
wires joined together with Y" or T junctions^^^^. 

In the present work, we examine a wider set of net- 
works constructed from one-dimensional wires and their 
junctions. The wires and the networks need not be con- 
fined to a plane. The wires may be bent out of the 
plane and may overpass each other, forming a three- 
dimensional network that is topologically distinct from 
a two-dimensional system. Fig. 1 illustrates a bilayer 
network of this type, though there is no barrier to form- 
ing an infinite stack of layers which is genuinely three- 
dimensional. The network is fixed in space and with 
permanent junctions, but we imagine that the electron 
densities in different parts of the network can be varied 
on a fine spatial scale by means of external gates, so that 
segments can be tuned at will between regions of topolog- 
ical superconductor (TS) and regions which are depleted 
of carriers and are, consequently, not in a topological su- 
perconducting state. For alternative ways to manipulate 
Major ana fermions in quantum wires, see Ref. 27 and 
28. Another possibility, which is conceptually interest- 
ing but may be more difficult to realize in practice, is 
suggested in the bottom half of Fig. 1. Here we envision 
a collection of flexible TS wire segments, which may be 
bent or tilted out of the plane at will, and which can 
be joined by junctions and redivided in different ways. 
For example, we imagine a process, illustrated in the fig- 
ure, where one end of wire 1 is attached to the middle 
of wire 2, and then the joint is broken in a way that 
one half of wire 2 remains attached to the end of wire 
1, while the other half becomes a disconnected segment. 
In either case, we require that the various wire segments 
are all proximity-coupled to a single bulk superconduc- 
tor, so that any differences in the phase of the induced 
superconducting order parameter between various parts 
of the network are strictly controlled. Further we require 
that the manipulations are done in such a way that the 
topological superconducting gap does not close. 

One of the unique features of a network of one- 
dimensional wires is that for an isolated wire segment, 
which is tunnel-coupled to the bulk superconductor but 
otherwise separated from the rest of the system, the 
parity of the electron number in the segment is con- 
served during adiabatic manipulations, and within the 
low-energy Hilbert space, the electron number parity can 
be related to the product of Majorana operators at the 
two ends of the segment. (The precise relation depends 
on sign-conventions, which must be carefully defined, as 
we discuss below.) If two occupied wire segments are 
adiabatically connected through a Y junction, and sub- 
sequently redivided, the total number parity will be con- 
served, but the parities of the individual wire segments 
may change. We show how the relations between elec- 
tron parity and Majorana operators can be exploited to 
predict the results of manipulations in which a Majorana 
site is moved through a Y junction connecting several 



wires. 

We also consider a variety of other manipulations 
that can alter the state in the zero-energy Hilbert space 
spanned by the Majorana operators. For the networks 
of semiconductor quantum wires discussed here, there 
are possible ways to manipulate the Majorana fermion 
state that do not have direct counterparts in a continuum 
two-dimensional system. In addition to the possibility of 
moving the Majorana locations in three dimensions, we 
find that the state may be manipulated by a rotation 
of the Rashba (spin-orbit-controlling) electric field (to- 
gether with the magnetic field) around an axis parallel 
to the wire, or a rotation of the Zeeman magnetic field, 
provided that the Zeeman field never points along the 
axis perpendicular to the plane containing the Rashba 
field and the wire direction. (Technically, the magnetic 
field should not come too close to this axis if one is to 
preserve the gap.) 

As in a two-dimensional network, the Majorana state 
may also be manipulated by means of controlled rotations 
of the superconducting phase difference between differ- 
ent parts of the network, which may be implemented by 
moving vortices through the bulk superconductor in re- 
gions between the superconducting nanowire segments. 
We also consider processes where field orientations are 
twisted along the length of a wire segment, or the wire 
itself is tied in a knot. 

We note that the behavior of Major anas passed 
through a Y junction and the behavior under rotations of 
the phase of the superconducting order parameter have 
previously been studied in Ref.^^ mainly using Kitaev's 
tight binding toy model^^. The latter model is topo- 
logically equivalent to the semiconductor-superconductor 
wire and a direct mapping exists for very strong Zeeman 
field^^. Exchange processes have been studied in details 
by Clarke, Sau and Tewari^^'^^. Our findings are in agree- 
ment with the previous works, although our methods are 
somewhat different. Our detailed study of the combina- 
tion of spin-orbit interaction with the Zeeman field allow 
us to analyze new types of manipulations, that are not 
possible in the context of Kitaev's chain model. Fur- 
thermore, our detailed study of the relation between the 
Majorana modes and the parity of the number of elec- 
trons within each wire leads us to conclusions regarding 
the physical observables associated with a series of braid- 
ings. The possibility of manipulating Majoranas in three 
dimensions has previously been discussed in a general 
way by Teo and Kane^^ and by Freedman et al.^^'^^, but 
not in the specific context of a wire network. 

We imagine initializing the system in a state where 
each wire has an even number of particles. This may be 
carried out if the wires are populated from full depletion, 
by Cooper pair tunneling from a superconductor. As 
electrons will enter the wire only in pairs, the resulting 
occupied section will be in a state of even parity. Fol- 
lowing the initialization we imagine a series of adiabatic 
manipulations, at the end of which the parity of the num- 
ber of particles in each wire is measured. Evidently, the 
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parity state does not depend on the exact location of the 
Majoranas at the end points of the wire and hence the 
manipulations are not required to correspond to closed 
trajectories of any kind. We find that the set of states 
to which the system may be brought by topological ma- 
nipulations is rather limited. After such manipulations, 
each wire may end up either in a pure state of even or 
odd parity, or in a mixture of even and odd parities, at 
equal weights. No other states may be constructed in a 
topological manner. 

Our analysis involves also a technical part. We show 
that the evolution of the system may be tracked by two 
equivalent methods. One method is based on a single- 
valued parametrization of the Majorana operators in 
terms of the time-dependent parameters of the system's 
Hamiltonian. We find that a single-valued parametriza- 
tion of the Majoranas necessarily involves discontinuities 
at arbitrarily chosen cuts. While these discontinuities are 
obviously absent from the physically observable quanti- 
ties, they are present in the relation between the par- 
ity operators and the Majorana operators. The second 
method uses a parametrization of the Majorana opera- 
tors that is multi- valued with respect to the parameters 
of the Hamiltonian. In this parametrization there are no 
discontinuities, but its being multi- valued makes it de- 
pendent on the history of the system. We show that the 
two methods are equivalent, but most of our analysis is 
carried out using the first one. 

The structure of the article is as follows. In Section II, 
we define the models to be considered, discuss the rela- 
tion of Majorana operators to the creation and annihila- 
tion operators for electrons, and define the sign conven- 
tions to be used in the paper. Specifically, in Subsection 
II A we discuss the general properties of the Majoranas 
in a wire with N spin resolved electron modes below the 
Fermi energy. In Subsection II B we consider the = 1 
case of a spin -polarized electron system coupled to a p- 
wave superconductor, and we discuss explicit solutions 
for the Majoranas at the end of a wire in that case. In 
Subsection II C, we analyze a quantum wire of electrons 
subjected to spin-orbit coupling and a Zeeman field, cou- 
pled by proximity to an 5- wave superconductor, which 
corresponds to A^ = 2 by our conventions. At strong 
magnetic field a projection of the N = 2 wire onto the 
low energy band reduces to the A" = 1 model. 

In Section III we explain the methods with which we 
analyze manipulations involving Majorana modes at the 
ends of a single wire. We first consider the general prin- 
ciples that determine the adiabatic time evolution of Ma- 
jorana fermion operators. In particular we show that any 
observable which is a product of zero energy Majorana 
operators is independent of time except at isolated in- 
stants when it may be multiplied by — 1 . We then analyze 
carefully the single- valued vs multi- valued representation 
of the Majorana end modes, and relate the Majorana op- 
erators to physical observables, in particular the parity 
of the number of electrons in each wire. Subsection HID 
gives examples for manipulations of one wire and their 




FIG. 1. (Color online.) Schemes for manipulating Majorana 
sites. Panel (i) shows two layers of a three-dimensional net- 
work of wires connected by Y junctions. Solid lines indi- 
cate wires in the first layer, dashed lines indicate wires in the 
second layer, and dotted lines indicate wires connecting the 
layers. All wires are tunnel-coupled to a bulk s-wave super- 
conductor with a definite superconducting phase. Gates (not 
shown) allow one to control the electron density in individual 
wire segments so that they can be either in a topological su- 
perconducting (TS) state or are depleted of electrons. Panel 
(ii) shows a portion of one layer with depleted regions (thin 
lines) and occupied regions in the TS state (thick lines). La- 
bels a-f show Majorana locations, which occur at the end 
of a TS segment or at a junction where three TS segments 
are joined. All wire segments must be very long compared 
to the coherence length in the superconducting segments in 
order for the interactions between Majorana variables to be 
neglected. Panels (iii)- (v) show an alternate scheme, where 
flexible TS wires may be joined together with a Y junction, 
and subsequently redivided in a different way. Labels a-d in- 
dicate Majorana positions. Panel (iii) shows an initial state 
where a is on the same wire as b, while c is joined to d. Panel 
(v) shows a final state where a is joined to d and b is joined 
to c. 

resulting transformations. 

In Section IV we discuss explicitly the possible sign 
changes that result when a Majorana is transferred 
through a Y junction, from one wire to another, as well 
as their physical significance. In Section V we use these 
results to study the sign changes and the physical con- 
sequences that occur when the locations of two Majo- 
ranas are exchanged by manipulation through a Y junc- 
tion. We examine both the case where the two Majoranas 
belong to different topological superconductor segments 
and the case where they are at opposite ends of a sin- 
gle topological segment. Results are summarized in Sec- 
tion VI. 



II. EXPLICIT MODELS OF TOPOLOGICAL 
SUPERCONDUCTING WIRES 

In this section we will review some explicit models for 
A'-mode wires that can form a TS phase. We will discuss 
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the structure of the zero energy state Majorana solutions 
at the end of the wires in several examples. 



A. N-mode quantum wires proximity-coupled to a 
bulk superconductor 

Majorana fermion excitations are known to appear in 
quantum wires which are proximity coupled to supercon- 
ductors^^^^'^^. Such systems are generally described by 
pairing Hamiltonians, 



H 



dw 



¥„{w)Ac,f,^Uw) + h.c. 



(1) 



These Hamiltonians are quadratic in the electronic cre- 
ation and annihilation operators and ^q,(i^), 
where a = 1, . . . , TV is a spin and/or channel index and 
w denotes the coordinate along the wire. (We will have 
occasion to consider wire orientations along arbitrary di- 
rections in space and reserve the notation x, z for the 
coordinates along the fixed laboratory coordinate axes.) 
The wire in the normal state is described by the single- 
particle Hamiltonian h^^ which includes the kinetic term, 
a scalar potential, as well as the spin-orbit and Zeeman 
couplings. The superconducting proximity effect is ac- 
counted for by the pair potential Aais which may de- 
pend on the momentum p^. While the single-particle 
Hamiltonian ho is hermitian, = ho, the Pauli prin- 
ciple demands that the pair potential is antisymmetric, 
A = -A^. 

The excitation spectrum of the Hamiltonian (1) as well 
as the associated quasiparticle (Bogolyubov) operators 
can be obtained with the ansatz 



dwK{w)^aH + <H^lH] (2) 



where the index a is summed over and the condition that 
jE be an eigenoperator of 7^;] = —Eje- This 

yields the Bogolyubov-de-Gennes equation Hi/j = Ei/j 
with the 2N x 2N Bogolyubov-de-Gennes Hamiltonian 



n 



ho 
At 



-h^ 



(3) 



and i/j = [u{w),v{w)]. Here, we collected the components 
Ua and Va into spinors u and v. 



As usual, the spectrum of the Bogolyubov-de-Gennes 
Hamiltonian is particle-hole symmetric. Indeed, 1-L obeys 
the relation H = —TxH^Tx, where Tx denotes a Pauli ma- 
trix in particle- hole space. Thus, if |?/;) is an eigenvector 
with energy then Krxlt/j) is an eigenstate with en- 
ergy {—E), where K denotes the complex-conjugation 
operator. Under certain conditions, the Bogolyubov- 
de-Gennes (BdG) equations possess zero-energy solu- 
tions, which correspond to Majorana fermions. These 
solutions satisfy (with an appropriate choice of phase) 
\ip) = Krx\^|J), which implies that u{w) = v^{w). In- 
troducing the notation ^^(it;) = u*^{w) for the Majorana 
solutions, the associated Bogolyubov operator becomes 



7 



j dw [g^{w)^^{w)^gl{w)^i{w)\ , (4) 



which satisfies the Major ana- fermion property 7 = 7^. In 
principle, a solution for 7 could be multiplied by an arbi- 
trary phase factor and still commute with the Hamilto- 
nian. However, the convention that a Majorana operator 
should be self-adjoint limits the allowable phase factors 
to ±1. 

The Hamiltonian H depends on a set of parameters 
such as the wire orientation or the direction of the ap- 
plied magnetic and electric fields. We will denote these 
parameters collectively as 11. In this paper, we will focus 
on processes in which these parameters vary slowly in 
time, denoting the initial parameters at time t = as fio 
and the final parameters at time t = T as ^It- We will 
analyze the resulting adiabatic time evolution of the Ma- 
jorana states given by g{w] fi) and of the corresponding 
Majorana operators 7(^2), focusing on two specific mod- 
els which we describe in the remainder of this section. 



B. Spinless fermions 

We first consider the Hamiltonian (1) for = 1, corre- 
sponding to spin-polarized electrons, with a pair poten- 
tial A = —ivpu)e^^^^^^ which originates from proximity 
coupling to a bulk px + Wy superconductor. Here, de- 
notes the overall phase of the order parameter A in the 
bulk superconductor, and a is the angle between the di- 
rection of the wire (given by the unit vector w) with the 
X-axis. (We assume for now that w lies in the x — y 
plane, so cos a = w • x, sin a = w • y.) We shall choose a 
branch cut such that — tt < a < tt. The Hamiltonian of 
the normal state includes only a kinetic term which we 
model with a quadratic dispersion so that 



H 



I 



dw 



^ ' ^ 2m 



*H + <b''{w)v{-ip^)e'^^+°'^<b''{w) + h.c. 



(5) 



This model of a one-dimensional spinless p-wave super- taev's toy model, with Bogolyubov-de-Gennes Hamilto- 
conductor is referred to as the continuum version of Ki- 
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man 



n = 



2m 



.e^^^+"V+ + /i.c.; (6) 



The matrices Ta denote Pauli matrices, which operate in 
particle-hole space. 

It turns out^^ that the continuum Kitaev model (5) 
also arises in the context of semiconductor wires with 
strong spin-orbit coupling in proximity to conventional 5- 
wave superconductors. The latter system, which we dis- 
cuss in the subsection II C, reduces to the continuum Ki- 
taev model in the limit of a strong Zeeman field. The ex- 
plicit connection between the models is given in Eq. (14). 

Majorana modes occur at the edges of the wire when 
the chemical potential within the wire is above the bot- 
tom of the electronic band (/i > in the case of quadratic 
dispersion). Substituting ip in Eq. (6) we find that 

the Majorana mode obeys the differential equation (here 
and throughout the rest of the paper h=l) \ 

-^9"{w) - ii9{w) - «e'(^+")5*'(«;) = 0, (7) 
Zm 

with the boundary condition ^(k; = 0) = 0; (we assume 
here that the wire ends abruptly near the origin of the 
coordinate system, and the occupied region has > 0). 
For hp > mv, the solution takes the form: 

7 = y dwg(w)^(w) + h.c] 
g{w) = e-^(["+^]^)/22v^4 



X sm 




(8) 



where /i = kp/2m, and [x]r is equivalent to x modulo 27r. 
This leaves an overall sign ambiguity in the definition of 
7, which we resolve, in the following, by choosing [x]r to 
lie in the interval 



-TT < [x]r < TT. 



Equivalently, we may write 



[x]r = [{x + tt) mod 27r] — tt. 



(9) 



(10) 



where "mod" is defined to give a value in the interval 
(0, 27r]. With these conventions, we see that for 6 = a = 

0, the first peak due to the term sin w^k'^ — (mvY is 

positive. 

Majorana states can also be localized in the bulk of the 
wire by allowing the local chemical potential /i to drop 
below the bottom of the band. Notice that since the 
chemical potential /i is defined as the difference between 
an electrochemical potential /io which is independent of 
position, and an electrostatic potential V{w) that may 
depend on position, spatial variations in /i are produced 
by variations m. V . At the interface where /i changes 



sign, a Majorana state will be localized. Now, we have 
to solve Eq. (7) separately for the region with negative 
chemical potential /iy < and the topological region 
with /iT > and match the wave function g and its 
derivative g' at the boundary between them. We assume 
that a + 6> = and choose units of mv'^ /2 for energy 
and mv for momentum. We assume that the topological 
phase, with jiT > 0, occurs for > and the non- 
topological phase (vacuum), with /iy < 0, for < 0. In 
these dimensionless units the equation for the Majorana 
zero modes becomes: 



-g'\w) - ^{w)g{w) - 2g\w) = 0, 



(11) 



with boundary conditions ^(0+) = ^(0~) and ^'(0+) = 
g'{0~). The full expression for the Majorana state, with 
a sign choice similar to the one defined above for the 
abrupt wire end, is 

g{w) = M~'^e~'^ ( cos ktw + — sin hct^o ) Q(w) 

\ KT J 

+ (12) 



Here we assume jir > ^ and /iy < and define A/" ^ = 

1 - /iy, 



1 > 0, 



and the unit step- function Q{w). If a + 6> 7^ 0, the ex- 
pression for g should be multiplied by e-*["+^]r/2 ^j^^ 
Majorana wave functions have a simple interpretation: 
in the "vacuum" the state decays with inverse decay 

length (y^I — 2iiy /mv'^ — l)mv l^^l-^^'^ y ^—2m/j.v 
In the topological phase the state decays with a de- 
cay length 1 / (mv) set by the induced superconduct- 
ing gap; on top of this, there are oscillations on the 
scale of the Fermi momentum in the superconductor 

tZT = mv^2iiT / ij^^'^) — 1 ^^^^^ > \/2m/j.T' Typical 
states at the end of the wire are depicted in Fig. 2. 

If initially the wire is along the x-axis then we have two 
Majoranas - one with a = and the other with a = tt. 
When we rotate the wire by 360° the term e^^t^+^l'^/^ 
will change sign for both Majoranas. However, a rota- 
tion of the wire by 180° will cause only one of them to 
cross the cut and hence only one of the Majoranas will be 
multiplied by a minus sign. We elaborate on this point 
is subsection HI C. 



C. Quantum wires with strong spin-orbit coupling 

Consider now the case N = 2^ corresponding to spin- 
ful electrons in the conduction band of a semiconductor, 
where the index a enumerates the spin orientations t 
and |. The electrons are subject to a Rashba spin-orbit 
coupling of strength originating from an internal elec- 
tric field that points in direction e perpendicular to the 
direction of the wire, a Zeeman field B, and a proximity- 
induced 5- wave pair potential A^^ = |A|e*^. 
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= -CX),/iT 
fly = -00,/Xt 



wmv 




wmv 



FIG. 2. Typical Majorana state at the end of the wire in 
the N=l case. In the upper panel the topological state is 
terminated in a vacuum with infinite potential while in the 
lower panel it has finite potential /Jv • The typical length scale 
for the oscillations and decay rate are depicted in the limit 
/XT ^ mv^ and /jy ^ —mv^. See Eq. (6) for the definition of 



In order to write the associated Bogolyubov-de- 
Gennes Hamiltonian, we find it useful to follow a mod- 
ified ordering of the Bogolyubov-de-Gennes spinor^^, 
namely tlj{w) = [u^{w)^ui{w)^vi{w)^ —v^{w)]. With 
this ordering, the off-diagonal particle-hole blocks of the 
Bogolyubov-de-Gennes Hamiltonian become diagonal in 
spin-space, 



n 



2m 



up ■ cr X e — /jj 



r^-B-o-+|A|e*V++h.c.. 



(13) 

Here, denotes Pauli matrices in spin space. In 
this representation, the particle-hole symmetry of 
the Bogolyubov-de-Gennes Hamiltonian follows from 
{H^—iTyT} = 0, where T = icjyK denotes the time- 
reversal operator and, with an appropriate choice of 
phase so that 7 is self-adjoint, a Majorana state will sat- 
isfy \^) = -iryT\ip). 

Of the parameters ft contained in this Hamiltonian 
that can be exploited to implement non-trivial transfor- 
mations within the ground-state manifold, our discussion 
will focus on the phase 6 of the superconducting pair po- 
tential, the direction w of the wire, the direction of the 
Zeeman field, and the unit vector e pointing along the 
Rashba field. 

A necessary condition for the wire to exist in a TS 
state is that the Zeeman field B has a non-zero projection 



onto the plane formed by e and w. Thus, we may always 
define a unit vector b pointing in the direction of this 
projection. For simplicity, we shall assume in most of 
our discussions that the component of B perpendicular to 
the (e, w) plane is zero, so that B = Bh. In the general 
case, a non-zero perpendicular component of B may be 
treated as an additional parameter of the Hamiltonian, 
along with, e.^., the magnitude of B and the profile of 
the variation in the chemical potential /i, which will affect 
the precise shape of the wave function associated with the 
Majorana operator, but will not be important for the sign 
conventions and electron parity manipulations we focus 
on in this paper. 

The Hamiltonian (13) can be reduced to simpler mod- 
els in a variety of limits^^'^^. In particular, it reduces 
to the Kitaev model discussed in the previous subsection 
in the limit of large Zeeman splitting. In this limit, the 
spin direction at the Fermi points is dominated by the 
Zeeman field, with only a small correction due to the 
spin-orbit coupHng. The strength of the s-wave proxim- 
ity effect is directly controlled by these small deviations, 
effectively causing p-wave-like correlations in the semi- 
conductor wire. At large magnetic field B mu^^ before 
the coupling to the superconductor, the two bands of the 
wire are well separated and only one band crosses the 
Fermi energy. If we project onto this band, the coupling 
to the superconductor gives rise to a BdG Hamiltonian 
of the form^^: 



n 



2meff 



P + Meff 



(14) 



with /ieff = fJ^ -\- B, '^eff = uA/B, 1/m* = l/m(l — 
mv? / B) and cos a = w • x. The 2x2 matrix operates in 
the particle-hole space (described by r matrices). This 
Hamiltonian represents a spinless p-wave superconductor 
in one dimension. Interestingly the phase of the order pa- 
rameter depends on the direction a of the wire - this fact 
will be useful when we will discuss interchange of Majo- 
ranas in a F junction with the wires oriented in specified 
directions on the surface of the bulk superconductor. 

The N = 2 case of semiconductor wires with strong 
spin-orbit coupling, subject to a Zeeman field, allows for 
a richer set of nontrivial transformations. In addition 
to rotations and changes of the phase of the order pa- 
rameter, one can consider variations of the direction of 
the Zeeman field B as well as the spin-orbit field ue. In 
addition, we consider rotations of the wire in three di- 
mensions under the assumption that the wire remains 
proximity coupled to the bulk 5- wave superconductor at 
all times. 

Away from the strong field limit, the Majorana so- 
lution of the BdG equations for = 2 has four com- 
ponents, which we write diS ip = -^(ipp^ip^) where the 
normalized two-component spinors i/jp and i/jh satisfy 
tjjh = T^l^p. Let us consider a wave function for a Majo- 
rana state, labeled by /, at the end of a TS wire segment. 
Since tjjh is determined by we need only describe the 
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FIG. 3. (Color online) Orientation angles defined in the text. 
The direction of the wire w in the fixed laboratory coordi- 
nate system is described by the angles a and (5. The electric 
field direction e, which is required to be perpendicular to the 
wire direction w, is rotated away from the (z, w) plane by an 
angle 5, not shown. Finally, 77 is the angle between e and the 
unit vector b, which is the direction of the projection of the 
Zeeman field B onto the e, w plane. In the text, we examine 
particularly the example where w is in the x — y plane, and 
e = b = z, so that /S — 7t/2, and 77 = ^ = 0. 

behavior of the latter. 

Let us suppose that the parameters such as the varia- 
tion in the chemical potential at the end of the wire, the 
magnitude of the spin-orbit coupling, and the magnitude 
of the Zeeman field 5, are fixed, along with the magni- 
tude of the superconducting pairing field |A|. The form 
of the spinor i/j^ will then depend on the orientation of the 

wire, the orientation of e and the orientation of b. For 
definiteness, we choose the wire orientation w to point 
from the depleted region into the TS region. Further we 
may define the variable w to be the distance along the 
wire taking as the origin the point where the chemical 
potential passes through zero, with positive w represent- 
ing the TS side. Let tljpo{w) be the spinor wave function 
in the case where w = x and e = b = z, with real A > 0. 
The function i/jpo is uniquely determined, except for an 
overall sign, which we choose using a convention that the 
a =t component has a positive real part at the point w 
where the wave function has its maximum amplitude [cf. 
the explicit forms for the wave function in strong mag- 
netic field limit, which were discussed in Eqs. (12) and 
(14) and illustrated in Fig. 2]. 

For any other choice of orientations, we may write 

^l{w',n) = Ui^po{w) (15) 

where Ui is a 2 x 2 unitary matrix. For brevity we will 
omit the subscript /, which enumerates the Majorana op- 
erators, throughout the remainder of this section. 

Let us first consider the case where b is parallel to 
e. (By definition, we have e ± w.) The orientation of 
the pair (w,e) relative to the original axes (x, z), may 
be described by an SO (3) matrix Rij, whose matrix ele- 
ments are the Cartesian components of the triad formed 
by w, (e X w), and e. The matrix can be described by 
specifying angles /3 and a which are the polar and az- 



imuthal angles of w (in the x, y,z coordinate system), 
plus an additional rotation angle S to describe the ori- 
entation of e relative to the (z,w) plane (see Fig. 3). 
Alternatively, one can characterize as a rotation by a 
specified angle about a specified axis. In general, R will 
have one unit eigenvalue, with a real eigenvector d, and 
two other eigenvalues of the form e"^*^. One can always 
choose A to lie in the interval < A < tt, and choose 
the sign of d so that R describes a rotation about d by 
a positive angle A. We may then define a unique SU(2) 
matrix S corresponding to R, by 

S = e-'^^-""'^ . (16) 

Note that if the matrix R is allowed to vary continuously, 
the matrix S will jump discontinuously whenever the ro- 
tation angle A reaches tt, as d will change sign at that 
point. 

Now let us introduce the superconducting phase an- 
gle ^, A = |A|e*^ . Then the matrix IJj describing the 
Majorana wave function may be written 

IJ = e'^I^S. (17) 

Again, we are faced with a choice of overall sign, which 
we must choose by some convention. We shall do this by 
defining the angle d to be in the range — tt < < tt. 

As an example, suppose that e is the z-direction. Then 
we must have /3 = 7r/2 and w = (cos a, sin a, 0). Then 

S = (18) 

Clearly, U will have a jump in sign if a crosses the cut 
at a = ±7r. Similarly it will have a jump if crosses the 
value ±7r. We may contrast this with the sign convention 
we employed for the = 1 wire. In that case, we had 
only a single cut, so the jump in sign occurs when 
reaches ±7r. 

Let us now consider the situation where b is not aligned 
with e. Since b is defined to lie in the plane of e and w, 
we may define an angle rj by 

e^^ = b • (e + zw) . (19) 

If we again assume, for simplicity, that w is oriented in 
the X — y plane, at an angle a relative to the x-axis, and 
that e is parallel to z, we can rewrite the Hamiltonian in 
Eq. (13) as: 

ho = — / - /i + ixpe-^^^"/V^e^^^"/2 
2m 

+ 5e"^^^"/^e"^^^^/V^e^^^^/^e^^^"/^ (20) 

This will clearly lead to an additional rotation of the 
Majorana spinor. If we generalize these formulas to the 
case where w and e are reoriented by an arbitrary SO (3) 
rotation, we find 

U = e^^/^S'e-^^^^/^ (21) 
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In order to make U single- valued, we shall require — tt < 
T] < 7T . This gives us an additional cut, with a jump in 
sign when r] crosses the value ±7r. 

In all cases, if the parameters of the Hamiltonian un- 
dergo a continuous change which returns to its starting 
value, so that fir = ^^o, then the matrix Uj will also 
return to its starting value. In doing so, it may have un- 
dergone a discontinuous sign change one or more times. 
The number Ng of these sign changes, mod 2, will be 
equal to the sum of the number of times that e*^ and A 
wind around the origin, plus one if the accumulated ro- 
tation of the pair (e, w) is equivalent to a rotation of 27r. 
We note that the net sign change, —1^% is a topological 
invariant of the path in parameter-space. It cannot be 
changed by any continuous deformation of that path so 
long as the gap remains finite throughout. This net sign 
change will be important in our analysis of the physical 
results of any manipulations of the Hamiltonian param- 
eters. 

To conclude, in this section we explored several topo- 
logical manipulations on the two Major anas at the end 
of a topological segment. Some of these manipulations 
do not involve any motion of the Major anas themselves. 
We showed that a rotation of the triad e,w,e x w by 
360° degrees around any axis leads to a multiplication of 
the Majorana operators by -1, provided the direction of 
the Zeeman field is simultaneously rotated so that it is 
fixed with respect to triad. As a particular case, if the 
magnetic field is kept parallel to the wire, and the wire 
is rotating around its own axis, so that e is rotated by 
27T while w and b are fixed, then the two Majoranas are 
multiplied by —1. Similarly, a change of the bulk s-wave 
superconducting phase, for example by taking a vortex 
around the wire, will lead to a multiplication of the Ma- 
joranas by -1. Additional factors of —1 may occur if the 
orientation of the Zeeman field is rotated relative to the 
triad, in such a way that the projection of the Zeeman 
field on the e — w plane winds around the origin. (This 
projection can never be too small or one would lose the 
gap in the wire). Majoranas can also be manipulated by 
depleting sections of a wire and pushing the Majoranas 
around bends. In section IV, we shall see how Majo- 
ranas may also be manipulated by pushing them through 
Y-j unctions. 



III. ADIABATIC MANIPULATION OF 
MAJORANA OPERATORS 

A. General principles 

In the previous section we introduced the zero energy 
Majorana states (or modes) at the two ends of each wire 
in the two models we consider. In this section we ex- 
plain the methods with which we analyze manipulations 
involving these modes. We first consider the general prin- 
ciples that determine the adiabatic time evolution of Ma- 
jorana fermion operators. Then, we relate the Majorana 



operators to physical observables, and emphasize their 
relation to the parity of the number of electrons in each 
wire. Subsequently, we apply the general methodology 
to manipulations of = 1 and N = 2 wires. 

Note that when analyzing the time evolution of the sys- 
tem under a time dependent Vt we will not limit ourselves 
to periodic trajectories. In a typical physical realization 
the requirement of a precisely periodic variation of the 
parameters, i.e., ^It = r^o, is rather stringent. When a 
vortex encircles a wire, or when a wire is rotated around 
its center, it is unphysical to require that all atoms in- 
volved in the motion return precisely to their original 
positions. Rather, we shall make the much less restric- 
tive assumption that the variation of 11 will be such that 
the Majoranas are always very far from each other, and 
the number of zero energy Majorana modes in the system 
is never changed. Under this assumption the Majorana 
operators 7^ keep their identity throughout their motion, 
and it is meaningful to track the time evolution of expec- 
tation values of products of Majorana operators such as 

{lilj) or {liljlkli)- 

According to Ehrenfest's theorem, if \^{t)) is the 
many-body wave function at time t, and A is an operator 
with an explicit dependance on the parameters 11, then 



d{^A\^) 
Jt 



dA 



~dt 



(22) 



If A is a product of zero energy Majorana operators cor- 
responding to the parameters fi, then [A, i^] = 0, and 
the first term is absent. Remarkably, the second term is 
also zero when ^ is in the ground state manifold, except 

^\ dis- 



A 



at isolated points in the trajectory where 

continuously changes its sign, because of our sign conven- 
tion. Other than these isolated points, the expectation 
value of A is time independent. This statement does not 
assume that the Majoranas involved belong to the same 
wire segment or to a single cluster of occupied segments 
connected by Y junctions, as long as they are connected 
to the same bulk superconductor, so that the supercon- 
ducting phase is well defined. 

To see that this is true, we first note that within the 
BdG model, we may write djj/dft = ^ic{ft)iji + Q, 
where Q where Q is a linear superposition of finite-energy 
Bogolyubov quasiparticle operators. If is located spa- 
tially far from all the other Majorana operators, so that 
their BdG spinors do not overlap, then the coefficients 
c{fl)i must vanish for all I ^ j. Moreover, if the param- 
eters are not at a point where the wave function for 
7j jumps discontinuously, then {^y j ^ j / d^l} = 0, since 
7I = 1. Thus we also have c{yt)j = 0. (See Supplemen- 
tary Material of Ref. [24] for more details.) On the other 
hand, if 11 crosses a special point lli where the BdG wave 
function for 7^ has a discontinuity in sign, then 7^ will 
also undergo a discontinuity in sign, without affecting the 
relation 7^ = 1. (Obviously, the derivative of 7^ is not 
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well defined at such a point.) Combining these results, 



we see that the expectation value ( ^ 



^ > is zero at 



dA 

all points except the points where the operator A changes 
sign due to a sign change of one of Majorana operators 
out of which it is composed. Notice that this conclusion 
is valid beyond the quadratic model of the BdG equa- 
tion, as long as the ground state subspace is spanned by 
hermitian Majorana operators 7^ that do not spatially 
overlap with one another. 

This observation leads to the conclusion that the set 
of parity states to which the system may be brought by 
topological manipulation is rather limited. To under- 
stand why it is so, assume that at the beginning there is 
a completely depleted network of wires. Then an appli- 
cation of gate potentials forms topological superconduct- 
ing sections with Majoranas at the end points. Assuming 
that during the application of the gate potentials there 
was no couphng between the sections we conclude that 
the parity of each section is even. 

Performing now topological manipulations such as 
braiding of Majoranas we can follow the evolution of the 
operator A = 7i7j. As argued earlier the only possi- 
ble operation is a multiplication of the operator observ- 
able (^4) = (^l^ll^) by a minus sign at the discontin- 
uous points. For example if 71 and 72 initially share a 
common section then they are in a pure even state with 
^ (7172) = 1. A braiding operation of 71 with 72 keeps 
the parity 1. Double braiding of 72 with an additional 
Majorana 73 will change the parity to —1. A single braid- 
ing of 72 with 73 leaves the wire section with 71 and 73 
at its ends. The initial value (7173) = acquires a mi- 
nus sign but remains zero. Hence the wire section must 
be in a mixture of either even or odd parities with equal 
weight. Generalizing this example for any pair of Majo- 
ranas we conclude that after the topological manipulation 
each wire section may end up in a pure state of even and 
odd parity, or in a mixture of even and odd parities, at 
equal weight. 

These observations point our attention to the im- 
portant, yet confusing, issue of sign discontinuities in 
the Majorana operators 7. We shall address this is- 
sue in detail, introducing two possible descriptions, one 
that avoids discontinuities at the price of multiple- 
valuedness, and one that uses discontinuities to ensure 
single- valuedness. In this discussion we will also re- 
late the expectation values of Majorana operators to the 
physically interesting observable - the parity of the num- 
ber of electrons in each wire. In an attempt for a clear 
exposition of the two descriptions of this problem, we 
start with a simpler, but related, problem: a spin 1/2 in 
a magnetic field. 



B. Spin in a time-dependent magnetic field 

We consider two quantum states of spin with s = 1/2 
in a magnetic field lying in the x — y plane, whose direc- 



tion forms a time-dependent angle a(t) with the x-axis. 
We analyze the phase the spin accumulates when a{t) 
winds adiabatically by 27r. The corresponding spinors 
have equal probability for spin cFz = ±1, which is analo- 
gous to the equal weight of the electron and the hole in 
Majorana fermions. The spin is described by the Hamil- 
tonian 



-B{cFx COS a ^ Gy sin a) 



(23) 



with B and a being parameters. The instantaneous 
eigenstates of i^cr depend only on the (fixed) a, and the 
ground state in the basis takes the form 



?s(a)) 



1 
71 



1 



(24) 



where x(<^) be chosen arbitrarily. When a varies 
slowly in time with initial angle ao, a spin that is initial- 
ized in the instantaneous ground state |gs(ao)) evolves 
in time into e*^^~*^^ |gs[a(t)]). The accumulated phase 
includes a dynamical part, Bt^ a Berry phase part. 



eB[a{t)]=lm 



/ ' 

J an 



da{gs{a) \dags{a)) , (25) 



and the contribution x [<^(^)] ~ X emerging from 

the explicit dependence of |gs[a(t)]) on x[<^(^)]- The dif- 
ference Xmon = xi^o + 27r) — xi^o) is, by definition, the 
monodromy phase. Notice that while the adiabatic evolu- 
tion of the ground state, including its accumulated phase, 
is uniquely determined by the Schrodinger equation, the 
way the geometric phase is split between the Berry phase 
and the monodromy depends on the choice of x(<^)- 

Indeed, a straightforward calculation of the Berry 
phase Ob{(^o + 27r), using Eqs. (24) and (25), yields 

Osiao + 27r) = TT + x(c^o + 27r) - x(c^o). (26) 

When computing the full geometric phase, Xmon — 
Ob {0^0 + 27r), the last two terms in (26) are canceled by 
the monodromy, resulting in a geometric phase of tt, in- 
dependent of the choice of x(<^)- 

Similar to gauge fixing, the choice of x(<^) is a matter 
of convenience, simplicity, and clarity of the calculation 
or argumentation. Three choices are worth discussing 
in the present context. If one chooses x(<^) satisfy 
x(q^ + 27r) = x{^)^ the instantaneous ground state is 
single-valued as a function of a and the total accumu- 
lated phase is given by the Berry phase Ob[o^o + 27r]. In 
contrast, choosing x(<^) = —a/2 makes the two com- 
ponents of the instantaneous ground-state spinors (24) 
complex conjugates of one another, not unlike spinor 
representations of Majorana fermion operators. For this 
multi- valued representation, the Berry phase vanishes, 
OB[a{t)] = 0, so that the geometric phase is accumulated 
entirely through the monodromy. 

Finally, it is also possible to construct a Majorana-like 
representation which is single-valued by choosing x(<^) = 
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FIG. 4. Plot of the function -[a]r/2, where [a]r = [{a + 
tt) mod 27t] — tt; cf. Eqs. (9) and (10). Choosing x(<^) — 

— [a]r/2 in Eq. (24) we have a single- valued Majorana-like 
representation of the ground state. 

— [a]r/2^ where [a]r is defined according to Equations 
(9) and (10); cf. Fig. 4. For this choice, x(<^) = x(<^ + 
27r), so that the total geometric phase originates from the 
Berry phase 6>b(<^o + 27r) only. However, the Berry phase 
is now accumulated in discontinuous jumps of x{^) ^it 
a = 7r(2n + 1) with n being an integer. We note that 
the time evolution of the physical wave function remains 
continuous, even for this discontinuous choice of x(<^) 



C. Application to Majorana variables 

As in this spin- 1/2 example, there is some freedom 
in defining the instantaneous eigenspinors tjj that cor- 
respond to the Majorana fermion operators, and there 
are analogs to both the single-valued and multi-valued 
choices we defined above. 

A zero-energy solution t/j of the BdG equation that 
starts in Majorana form (ipp^Ti^p)^ at t = stays in 
that form throughout its time evolution, and satisfies 
{ilj\dt^l^) = at all times. As a consequence, upon com- 
pletion of a periodic trajectory, = 11^, both the zero- 
energy Majorana spinor t/j and the corresponding Ma- 
jorana operator 7 must return to their initial values up 
to a possible sign. As before, if the Majorana operators 
are defined singly- valued as a function of 17, they acquire 
discontinuous minus signs at particular values of ft. If 
continuity is preferred and the price of multi-valuedness 
is paid, then for every value of ft the Majorana opera- 
tors are defined only up to a minus sign, and these minus 
signs depend on the entire history fl{t). 

1. Electron number parity. 

Let us now relate this choice of Majorana operators 
to physically measurable quantities, in particular to the 
electron parity. We start with the simplest case, that of a 
single wire, with Majorana modes 71,72 at its two ends. 
For such a wire we may define an electron parity operator 
P, which may take the values ±1, depending on whether 
the number of electrons in the wire segment is even or 
odd. We naturally expect that within the zero-energy 



Hilbert space, the parity operator should be related to 
^71 72, with the eigenvalues being ±1. Far less clear, and 
far more confusing, is the assignment of the two eigen- 
values to even and odd number of electrons, and the way 
this assignment evolves when Q, varies in time. 
Generally, we may write 

P = ikjij2 , (27) 

where /c = ±1. The value of k may depend on the ori- 
entations of the wire ends and the sign conventions used 
in the definition of the Majorana operators. It is easy 
to see that different choices for the dependence of 71 , 72 
on ft lead to different values for k. We note, however, 
that the actual electron number parity cannot change in 
manipulations of a single wire, as the parity is an in- 
variant quantity in any adiabatic manipulation that does 
not bring it into contact with another wire and does not 
deplete the wire to the point where all electrons are ex- 
pelled. In turn, this reflects the fact that the tunnel 
coupling to the bulk superconductor allows only electron 
pairs to tunnel between the superconductor and the wire, 
unless one supplies enough energy to inject an electron 
or hole above the energy gap of the bulk superconductor. 

Consider a straight wire lying in the x — y plane, at 
an angle (j) relative to the x-axis, with e and b in the 
z-direction. We define angles aj as the azimuthal angle 
relative to the x-axis of a line parallel to the wire at 
end j = 1,2, directed from the depleted region into the 
region occupied by the topological superconductor, with 
the restriction — tt < aj < tt. Then, ai = (j) and a2 = 
(j) ±7r. The Majorana operators 71, 72 depend on ai, 0^2 
respectively, and that dependence may be chosen to be 
either single- valued or multi- valued with respect to aj. 
We imagine now that the wire is rotated by an angle tt 
around its center, and examine what happens to P, /c, 71 
and 72. This is a point where the difference between 
the single- and multi-valued representations should be 
carefully followed. 

If the Majorana operators are defined in a way that is 
single- valued with respect to Ofj, then the interchange of 
ai with a2 leads to the interchange 71 72 and 72 ^71. 
That is, after the interchange, the expression for 71 in 
terms of the electron creation and annihilation operators 
is the same as the expression for 72 before the interchange 
and vice versa. Since ^72 71 = —^7172 and since the num- 
ber parity is not changed by the rotation, the pre-factor 
k has necessarily changed sign. That change of sign takes 
place at the value of (j) where one of the 7's acquires a 
discontinuous sign change. The parity operator as ex- 
pressed in terms of creation and annihilation operators 
for the electrons does not change discontinuously at this 
point. Similarly, the actual physical value of the parity 
does not change during the manipulation. 

In contrast, if the Majorana operators are defined in a 
way that is multi- valued with respect to a^, they gener- 
ally depend on e^^*"^/^. Then, a rotation by tt leads to 
the transformation 71 ±72 and 72 TTi? with the 
proper sign being determined by the sense of the rotation. 
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With that choice, the parity operator stays P = ^7172 
and k does not change. Again, no discontinuous change 
occurs in the parity operator when expressed in terms 
of electrons creation and annihilation operators or in the 
physical value of the parity. 

Comparing these two descriptions, the advantage of 
the single- valued choice of the 7's is that all quantities 
are fully determined by the instantaneous values of the 
parameters, in this case a^, and are independent of the 
values of the parameters in earlier times. The depen- 
dence on the history of the system - preparation, manip- 
ulation etc. - is all in the many-body wave function of 
the wire. These advantages come at the price of having 
discontinuities in the 7's and in the way P is expressed 
in terms of the 7's. These discontinuities appear in arbi- 
trarily chosen locations in parameter space. In contrast, 
the advantage of the multi- valued way is in the absence 
of any of these discontinuities, while the disadvantage is 
in an arbitrariness in the initial definitions of the oper- 
ators 7i , 72 and the need to follow the time evolution of 
these operators to be able to relate them to the parity P. 

The two descriptions may be extended to more compli- 
cated situations. Generally, any N = 1 wire may be ob- 
tained from an A/' = 1 straight wire by continuous bend- 
ing. In the single-valued description, as the wire is bent 
from being straight to its final form, sign discontinuities 
may occur in the Majorana operators 71,72- These dis- 
continuities occur when either of the angles aj crosses a 
cut. Whenever such a sign change occurs, a similar sign 
change occurs in the prefactor thus guaranteeing that 
the parity of the number of electrons remains continuous. 
It should be noted that the form of the parity operator 
is not fully determined by aj. Two wires with the same 
values of aj may have opposite values of /c, if the defor- 
mation of one into the other involves an odd number of 
crossings of cuts. 

In the multi-valued description, the prefactor k does 
not vary at all with the bending of the wire, but the 
values of the angles a^, of the factors e^^^l^ and hence of 
the operators 71, 72 depend on the trajectory that brings 
the wire from its initial to its final form. 

For the N — 2 case, two wires may be deformed into 
one another as well, but the deformation is more compli- 
cated, since the Majorana modes depend not only on the 
directions of the wire's ends, a^, but also on the direc- 
tions of e and b. Again, in the single-valued description 
a discontinuous sign change associated with cuts in the 
definition of 7^ as a function of 11 leads to a sign change 
in /c, leaving the parity continuous. And again, the value 
of k is not fully determined by the values of , e and b 
at the two ends. Rather, a variation of these parameters 
along the wire may lead to sign changes in k relative to its 
value for a reference straight wire where, say, = 0, tt 
and e, b are parallel to the z-axis. 



D. Examples 

We now analyze two A^ = 1 examples, of a straight and 
a curved wire. For both, we calculate the parity operator 
in a representation that is single-valued with respect to 
ri. Following that we will discuss in subsection HID 3 
the N = 2 case. 



1. Straight wire 

To find the parity operator for a straight wire we use 
the discrete version of the Hamiltonian of Eq. (5) with 
t = 11^ o? = l/(m/i), A = v/a and the number of sites 
A' = L/a-\-l. Here with a being the lattice constant and 
L the wire length, the Hamiltonian becomes 



H 



N-1 

1^1 



I H+i 



+ h.c, (28) 



where a describes the orientation of the wire in the x — y 
plane, directed from the end / = 1 to the end / = A". 
Since, in our conventions for Majorana operators at the 
end of a wire, we have defined the angle a to describe 
the orientation of a vector pointing into the wire from 
the end, we note that a = o; for a Majorana at the end 
/ = 1, but a = q; ± TT at the end / = A". 
Let us now define 



ai = e'^l^ci 



hi = i{e 
which gives 



ci 



-[e + a]r, (29) 



iaihi = 1 - 2c\ci 
{auam) = {buhm} = 28im, {bi,am} = Q- 



(30) 



If we now consider the case | A| =t > 0, the Hamiltonian 
may be rewritten as: 



N-1 



H = t^ iai^ibi 



(31) 



1^1 



The Hamiltonian clearly commutes with ai and with 
bri' These two operators are therefore the Majorana end 
modes, up to possible factors of ±1. In order to deter- 
mine the sign, we compare the definitions in (29) with the 
sign conventions defined for the continuum wave function 
(8). We see that the Majorana operator at / = 1 may be 
written as 7^ = ai, but at the other end, we have 

75 = 6ivsgn([^ + a],) . (32) 

In terms of the Majoranas, the parity operator is 

N N 



P = Yl{l-2c]ci)=]l{iaik). 



(33) 
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When the wire is empty ah sites are empty and the 
parity is 1. In the ground states of Hamiltonian (31) 
{ibiai-^i) = 1 so that for the ground state manifold the 
parity operator is 

N-l 

P = iaibN Yl = ^^i^iv = ikjalb^ (34) 

k = sgn{[0^a]r) . (35) 

Although we have explicitly considered only the special 
case where t = \ A\ and the chemical potential is at zero, 
the relation between P and ijalb cannot change discon- 
tinuously if the parameters are varied, as long as there is 
an energy gap in the bulk of the wire, and the orientation 
of the wire is fixed. Furthermore, when the chemical po- 
tential is not fixed at the center of the band, we can take 
the continuum limit of a ^ 0, with the electron density 
held fixed. This recovers the continuum = 1 model 
defined earlier. 



2. Curved wire 

To determine k explicitly for a curved wire in two 
dimensions we have to follow the configuration of the 
wire in the two dimensional space. Assuming that the 
curve describing the wire is {x{s)^y{s)) we define an an- 
gle a{s) as the orientation of the tangent vector at point 
s, choosing the vector to point along the wire direc- 
tion leading from a specified end {s = 0) to the other, 
and requiring that a is a continuous and slowly vary- 
ing (on the coherence length scale) function of s. We 
see that tana = dy/dx and the curvature n = da/dw = 
{x' y" — y' x") / ix''^ ^y''^) with w being the arc length coor- 
dinate along the wire, defined by dw = Y^^^^^"^r"^<^<^- We 
define the net orientation change as Aa = a{L) — a(0) = 
tzdw^ where L is the length of the wire segment and 
we have chosen the parameter s to be the arc length w. 
Note that \Aa\ can be larger than 27r if there are loops 
in the wire. We shall assume here that the phase of the 
bulk superconductor is a constant along the wire. 

If |Aa| lies in the interval 27r(M- 1) < |Aa| < 27r(M+ 
|), (with M being an integer), then our result for the 
factor k relating the electron number parity operator to 
the product ijalb is 

k = {-l)^q , g = sgn ([^ + a{L)]r - [0 ^ a(0) - n]r) . 

(36) 

To derive this result, consider what would happen if 
one were to adiabatically change the parameters in the 
Hamiltonian so that the wire segment becomes straight, 
while holding fixed the orientation a(0) at the first end. 
The wire orientation a{L) at the second end will change 
continuously from its initial value ao{L) to a final value 
equal to a(0). Of course the expectation value of the 
electron parity will be unchanged in the procedure. 



During the process of straightening, the value of M 
will change by ±1 whenever a{L) — a{0) = (2n + l)7r, 
but the sign of q will also change at these points. The 
sign of q will also change when (j)-\-a{L) = (2n + l)7r, but 
the operator 75 at the end w = L will also change sign 
at such points. Therefore the product ik{ja7b) will be 
unchanged during the process. When the wire is straight, 
however, we have Aa = 0, and we see that [0 + Q;(L)]^ 
and [6 + a(0) — 7r]r have opposite signs. Therefore (36) 
reduces to (35). 

3. Wires with N = 2 

For N = 2 the Major anas depend on a richer set of 
parameters, and thus cross more types of branch cuts. 
For example, cuts may occur if we consider situations 
where the direction of e is allowed to change along the 
wire, or if the wire orientation is lifted out of the plane, or 
if the angles rj or are allowed to vary along the wire. In 
these cases we must follow all the orientation parameters 
along the wire. 

To describe the general situation, we begin by defin- 
ing the unit vectors e, b, w, at all points along the wire, 
where w points in the wire direction from end a to end b. 
We require that e is perpendicular to w at all points, and 
we define r]{w)^ as before, as the angle between e and the 
projection of b onto the plane formed by e, w. We define 
an SO (3) matrix R{w) which describes the orientation 
of e,w relative to the laboratory axes x, y,z. We may 
define angles A77 and A^ to describe the net change in 
T] and 6 along the wire, when these angles are defined to 
vary continuously as a function of w. We may also define 
an SU(2) matrix S{w) as the matrix corresponding to 
the relative rotation R{w)R~^{0), with the requirement 
that S{w) is a continuous function of w. With the sign 
conventions for 7^ and 75 specified earlier, we find that 
the expression for k has the form 

fc = (-l)^gsgn[tr5(L)] , (37) 

M = Int(A7//(27r)) + Int(Al9/(27r)) , (38) 

where Int(x) denotes the closest integer to x, and Q = ±1 
depends on the matrices R{0) and R{L). Although we 
have not found a simple form for the function Q, the es- 
sential point is that its value may be computed if desired 
in any instance, and it depends only on the end points of 
the orientations and not on the path between them. 

To derive the result we may again begin by considering 
a straight wire in the x — y plane, with b independent 
of position and e = z. We again define the angle a 
as the orientation of the wire directed from the end with 
Major ana operator 7^ towards the end with 75. If we now 
assume that we are in the strong- field limit, we can map 
the problem onto the A^ = 1 Kitaev model considered 
previously. The constant k in this case is given by 

k = sgn{[a]r) . (39) 
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The difference between this expression and (35) arises 
from the different convention we have adopted for the 
branch cuts in the two cases. In the present N = 2 case, 
the jump in sign of a Majorana variable occurs when the 
orientation of a (pointing into the wire from the end) 
crosses tt mod 27r, independent of the phase of the super- 
conducting order parameter. The result is also indepen- 
dent of the angle r]. 

The results for a bent wire of arbitrary orientation, 
and with r] and that vary along the wire, may be ob- 
tained by deforming the wire adiabatically into the case 
of a straight wire in the x — y plane. The value of k will 
be multiplied by —1 each time the orientation parame- 
ters of the Majorana at either end of the wire cross a 
discontinuity surface. 

As a particular example, consider a bent wire in the 
X — y plane, with e = z everywhere. Suppose that a{w) 
describes the local orientation of the wire, as in the sub- 
section on TV = 1, and suppose the total bending angle 
Aa is less than 27r. In this case we find 

k = sgn {[a{0)]r - [a{L) - 7r]r) = sgn {[aa]r - [c^b]r) , 

(40) 

where aa and describe the orientations of the Majo- 
ranas at the two ends. This result coincides with the 
result (36) for the N = 1 case when the phase of the 
superconducting order parameter is zero. 



IV. MOVING MAJORANAS THROUGH A Y 
JUNCTION 

In the previous section we analyzed two Majorana op- 
erators at the ends of a single wire, and showed how their 
states may be manipulated in such a way that preserves 
their identities but implements non-trivial transforma- 
tions. In this section we enlarge the system to include 
several wire segments connected by junctions, which will 
be necessary in order to perform non-trivial transforma- 
tions on the state of the system, and in particular, braid- 
ing. 

Initially, our system consists of an even number of Ma- 
joranas, distributed into pairs, with each pair consisting 
of end modes of the same topological segment. Any non- 
trivial transformation will require a rearrangement of the 
pairing configurations while keeping the identities of the 
Majorana states distinct, so that the quantum informa- 
tion encoded in their state is preserved. Thus, the pro- 
totypical situation we will eventually consider is one in 
which at t = we have two topological segments, with 
the Majorana operators 7a, 76 located at the two ends of 
one segment and 7c,7(i located at the ends of the sec- 
ond segment. At the final time t = T, one segment has 
7a, 7c at its ends and the other has 75, 7d at its ends. 
This exchange of the two Majoranas should be carried 
out through an adiabatic evolution of the parameters 11, 
without introducing a coupling that will modify the en- 
ergy of any of the Majorana states. The analysis of such 



manipulations requires care, as exemplified by the fol- 
lowing scenario. If both wires are initially in a state of 
even parity, the interchange of 75 with 7c will result in 
a state where the parity of both wires has a zero ex- 
pectation value. If the interchange is then repeated, the 
two wires may end up both at even parity, or both at 
odd parity, depending on the precise details of the inter- 
change process. It is the analysis of this process which 
we now outline. We shall accomplish this process using 
a series of elementary manipulations. The elementary 
tool for these manipulations is the Y junction. At the Y 
junction three wires come together each of which may be 
"empty" or "full" {i.e., in a TS state). If the number of 
TS segments at the junction is odd, there will generically 
be one zero-energy Majorana state at the junction. If the 
number of wires is even, there will be no such state. 

The case of three full (i.e., topological) wires may be 
understood by depleting each wire just before it reaches 
the junction. If tunneling between the wire ends were 
completely forbidden, there then would be a zero-energy 
Majorana state at the end of each wire, i.e., three Majo- 
ranas all together. However, when tunneling is allowed, 
the three states will be split, with two combinations of 
the original Majorana operators pushed to finite energy, 
and one combination forming a new Majorana operator 
with precisely zero energy. 

To see this formally, we use the Hamiltonian de- 
scribing the tunneling between Majoranas 1,2, and 3 
at the ends of the three wires near the Y junction, 
H = -^"^Zi m^i 2 s^irnjnm- The real matrix elements 
tnm present the coupling between the Majoranas. Since 
H = we require that tim = —tmi- Writing the 
solution for the equation [H, Tn] = —EnTn as Tn = 
Vniji + Vn2j2 + VnSjs^ wc obtain the matrix equation 
for the amplitudes Vrnn- 




En 



■ (41) 



This equation ha s two particle-h ole symmetric solutions 
with energies i\/^i2 + ^23 + ^135 ^^^^ energy so- 

lution corresponding to the true Majorana located near 
the junction. The zero energy solution is Vorn ^ Qkmtik 
(with eikm the Levi-Civita symbol). Explicitly, its nor- 
malized form is 




(42) 



If the tunnel coupling between two of the wire ends 
is much stronger than the coupling to the third (e.g. 
ti2 ^ ^23, ^13), then, as can be seen explicitly by examin- 
ing the vector vq in this limit, the zero energy Majorana 
operator will be primarily located at this third wire end, 
while the first two wires combine to form, effectively, a 
single wire, with a finite energy fermion state inside the 
superconducting gap of the TS sections. 
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FIG. 5. (Color online) Moving a Majorana across a Y junc- 
tion, (i) Basic set up of a y junction. The bold (blue) line 
represents a topological superconductor (TS) sector while the 
thin (red) line denotes an empty wire. The three branches of 
the junction, 1, 2 and 3 are at angles < 01, 02, 03 < 27r re- 
spectively to the X axis. (For the N=l case these are also the 
phases of the effective p-wave superconductor). In this panel, 
a Majorana a is located on branch 2. The configuration of the 
wires near the junction is denoted by the symbol ^Y, where 
the bar symbol above the Latin letter a denotes that it is the 
end of a TS sector stretching away from the junction. Ma- 
jorana a is characterized by the angle aa = 02- (ii) A basic 
operation is to move Majorana a from branch 2 to branch 3. 
To do so we first move the Majorana to the Y junction, which 
is done by adjusting gate voltages to fill the remainder of wire 
2 with electrons in the TS state, (iii) Next we bring a portion 
of wire 3 near the junction into the TS state. This puts the 
Majorana on leg 3. The symbol t in the Va denotes that the 
entire branch 2 is topological and has no Majoranas. The bar 
has been removed from the letter a because the Majorana is 
now at the end of a TS segment connected to the Y junction. 
It is characterized by an orientation aa — 03 H-tt. Lower pan- 
els (i')-(iii'). Shifting a Majorana from branch 1 to branch 3 
when both branch 2 and 3 are in the TS state. We denote 
this by "y/^4V^. As shown in the text, this operation gives 
an additional sign change relative to the case in the upper 
panels. 

The basic manipulations at a F junction involve push- 
ing a Majorana state from one wire to another through 
the junction. Here we may distinguish two situations, 
illustrated by the upper and lower panels of Figure 5. In 
the upper panels we have one full wire coming in from far 
away (wire 2), while the other two are empty. Initially, 
the full wire stops short of the Y junction, and the Ma- 
jorana state at the end of the filled region, is located on 
wire 2. At the end of the manipulation, wire 2 is com- 
pletely full, and wire 3 is filled to some finite distance 
from the junction. The Majorana, at the end of the filled 
section, is now found inside wire 3. This manipulation, 
which does not involve a coupling between two separate 
topological segments, is rather simple. 

The lower panels describe a more complicated manip- 
ulation. In this manipulation we again transfer a Majo- 
rana from wire 2 to wire 3, but here we have the opposite 
situation, where all three wires are full far from the junc- 
tion. In the initial state, we assume there is a depleted 



region in wire 2 between a point r2 and the junction, 
while wires 1 and 3 are everywhere full. The Majorana is 
located at the point r2. At the end of the manipulation, 
the depleted region has been moved to wire 3, and the 
Majorana is located in wire 3. The manipulation involves 
a transfer of a Majorana between wires. 

As long as we change the parameters of our Hamil- 
tonian at a rate that is small compared to the energy 
of the lowest positive energy states that occur in the 
course of the manipulation, positive energy states will 
never be occupied, and we may ignore them in our con- 
siderations. For the analysis of the time evolution of the 
single zero energy Majorana operator throughout the ele- 
mentary manipulations, we find it convenient to revert to 
a representation in which the Majorana operator evolves 
continuously with the parameters of the Hamiltonian. 
We denote this representation of the Majorana operator 
by 7(^1). As discussed in Sec. Ill A, in this representa- 
tion, the expectation value of the product of 7(11) with 
any number of Majorana operators far from the junction 
will be independent of time. We can always choose the 
sign of 7(ri) to be consistent with our sign conventions 
at the beginning of the manipulation, (t = 0), when the 
Majorana is on wire 2. However, we are not guaranteed 
anymore that 7(11) is still consistent with the conven- 
tion at the end of the manipulation (t = T) when the 
Majorana is located on wire 3. If the conventions are 
not satisfied, we will encounter a factor of -1 when we 
convert 7(11) to a Majorana operator obeying our single 
valued convention, and these factors of -1 must be in- 
cluded when we compute expectation values of physical 
operators. 

Fortunately, the rules for when one needs a factor of 
-1 are independent of most details of the Y junction. In 
the remainder of this section, we show how the rules can 
be deduced from the requirement that the total electron 
number parity be conserved, together with our previously 
developed relations between the electron parity operator 
of an occupied segment and the product of Majorana 
operators at its ends. 

In the following analysis, we discuss the behavior of a 
network of semiconductor wires with N = 2. We assume 
that near the junction the phase of the superconductor 
and the angle rj are constant, and the wires are in the 
X — y plane, with e in the z-direction. (Results for any 
other orientation of the Y junction can be obtained by 
performing a rotation of the system and counting the 
number of times the parameter set ft at the location 
of a Majorana crosses a discontinuity surface during the 
rotation.) 



A. Y Junction Notation 

We first define a notation to keep track of the vari- 
ous arrangements of TS wire segments and Majoranas 
in the vicinity of a junction, such as those exhibited in 
the various panels of Fig. 5 . We mark Majorana states 
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on the branches of the Y with a lower case Latin letter 
located about a Y symbol, at the appropriate location. 
If the Majorana state is an 'in' state (it is at the end 
of a topological segment which stretches away from the 
intersection), we mark the state with an extra bar. The 
presence or absence of a bar leads to a difference of tt in 
the value of the angle a describing the Majorana orien- 
tation. When a wire branch is filled in the TS from the 
junction all the way to "infinity" , so there is no Majorana 
near to the junction on that branch, we indicate this by 
including the letter t on the corresponding branch of the 
Y symbol. 

For example, consider: 

^Y. (43) 

which describes the state in Fig. 5(i). Wires 1 and 3 are 
empty, while branch 2 contains the Majorana state 7a, 
terminating a topological segment that stretches away 
from the intersection. Another example is: 

'Y,. (44) 

which encodes a Majorana state on branch 3, which ter- 
minates a topological phase that stretches to infinity on 
branch 2, as the t indicates. This is the situation in 
Fig. 5(iii). 



B. Derivation from Number Parity Conservation 

Let us first consider the situation illustrated in the 
upper panels of Fig 5. We suppose wire 2 is straight and 
long but finite, so that it has an additional Majorana at 
its far end (not shown in the figure) which we label 72. 
The orientation of the Majorana 72 is a2 = 02 i^r, while 
the original orientation of 7^ is 0^2 = ^2 • 

In our notation, the Majorana operator keeps the same 
label a throughout the manipulation. However, the rela- 
tion of the operator to the fundamental electron creation 
and annihilation operators will be different before and 
after the manipulation, as the parameters 11 will be dif- 
ferent. To emphasize this difference, we shall use the 
symbol 7^, in this section, to denote the Majorana oper- 
ator in the final state in the single- valued convention, and 
use a'^ to denote the corresponding orientation. In the fi- 
nal state illustrated in panel (iii) of Fig 5, the orientation 
of Majorana a will be a'^ = (/)3 ± tt. 

In the initial and final states, the number parity oper- 
ators are respectively given, according to (40), by 

kij = sgn ([0i - 7r]r - [(/)j]r) . (46) 

We now relate the operator 7^ to the continuously 
evolving operator 7(^2). Since the expectation value 
(727(11))^ in the final state is equal to the expectation 



value (727a)o ^^e initial state, and since the expecta- 
tion value of the number parity must be the same in both 

states, /C23 {l2la)T = ^22 {^21(S^))t' ^^^^ ^ 

in the final state, the operator 7^ is equal to 8237(^2), 
with S23 = /c22^23- The product /c22^23 will equal -1 if 
and only if the x-axis lies within the area between two 
lines of orientation ^3 and (/)2 — tt, with the restriction 
that the opening angle have magnitude < tt. 

In accord with the discussion in Section III A, we know 
that if A is any operator constructed as a product of 
Majorana operators far from the Y junction in question, 
the expectation value {A^a{^)) will not depend on time. 
We therefore see that 

{Ai^l = S23(A7a> (47) 

where the expectation value of the left hand side of the 
equation is taken at time T, and that of the right hand 
side is taken at time zero. 

The inverse of the above manipulation begins with a 
TS that fills the entire wire 2 and extends, through the 
junction, part way into wire 3, so that there is a Majorana 
located in wire 3. If we now adiabatically deplete the 
occupied region of wire 3, along with a region close to 
the junction in wire 2, we will have moved the Majorana 
from wire 3 to wire 2. Clearly the sign change associated 
with the inverse process, which we denote 532, will be the 
same as ^23. 

When all three wires are in the TS state far from the 
junction, as in the lower panels of Fig. 5, the analysis is 
slightly more complicated. We assume that all wires are 
straight, and we include two more Majorana operators, 
71 and 73 at the end of their respective wires. In this case, 
the number parity for individual wires is not conserved, 
but the total parity is. The parity operators in the initial 
and final states are given by 

= -k22k3ll2lal3ll P^ = -^33^2l737a727l (48) 

Taking into account the anticommutation rules for Ma- 
jorana operators, we know that the expectation value of 
737(r2)727i in the final state should be equal to the ex- 
pectation value of —727a 7371 in the initial state. This 
tells us that for ft = I^t, being the parameters at the 
final state, the operator 7^ is equal to 5237(^2^), with 
-^23 = —^22^31^33^21- Furthermore, if A is any product 
of Majorana operators far from the junction, we have 
{A'y'^)T = 523(^7^)0, where the subscript indicates the 
time at which the expectation value is taken. 

As an example of the above, consider the case where 
(^1,^2,^3) = (7r/6,57r/6,37r/2). Then /C22 = -l,/^3i = 
1, kss = 1, A:2i = 1. In this case we have no sign change: 
in the final state we have 7^ being equal to the con- 
tinuously evolving 7 at time T. On the other hand, if 
we rotate the system by — 7r/3, so that (0i, (/>2, (/>3) = 
(37r/2,7r/6,57r/6), we find /C22 = l,/^3i = -1,^33 = 
— 1,^21 = 1, so in the final state we will have 7^ = 

These results can be readily generalized to the case 
of arbitrary angles between the wires. We find that the 
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factor of -1 occurs if and only if the negative x-axis lies 
within the area between two lines of orientation cj)^ and 
02, with the restriction that the opening not include the 
direction — tt. 

The results derived above can be stated in another 
way. For both cases considered in Fig. 5, the necessity 
for including a minus sign at the end of the manipula- 
tion can be determined by defining a continuous path for 
change in the orientation angle between the original 
orientation in wire 2 and the final orientation in wire 3. 
If this path crosses the cut at angle tt mod 27r, we must 
multiply the Majorana operator by a minus sign. In the 
case where one incoming wire is occupied far from the 
junction (upper panels of Fig. 5 ), we must define the 
path between the original and final angles by requiring 
that a is never equal to the wire orientation angle (j)2' 
If we think of the continuous path as a set of Majorana 
orientations which would occur if one continuously bent 
the end of wire 2 until it aligned with wire 3, the condi- 
tion effectively says that the wire would not double back 
on itself. In the case of three full wires, illustrated in the 
lower panels of the figure, the path should be chosen by 
requiring that a must pass through the orientation angle 
of the third wire at the junction. 

If the Y junction is not oriented parallel to the x- 
y plane, or if the orientation of b is not parallel to z, 
the factors kij in the parity operators will no longer be 
given simply by Eq. (46). Nevertheless, the relations 
= ^22^23 and ^23 = —^^22 ^31 ^33^21 Still hold, as they 
followed directly from conservation of overall parity. 

V. MAJORANA EXCHANGE 

The process of exchanging two Major anas can now 
be performed by compounding elementary processes of 
the type described above, where a Majorana is moved 
through a Y junction from one wire to another. 

An exchange process has two distinct fiavors: either 
the exchange is between two Majorana states at the edges 
of the same topological segment, or between states be- 
longing to different segments. A configurations of each 
type is illustrated in Fig. 6. 




(i) 3 (ii) 3 

FIG. 6. Two initial states for a Majorana exchange process, 
(i). Two Majoranas on the same TS segment, (ii). Majoranas 
on different TS segments. 

When the Majorana states belong to the same topo- 
logical segment, an exchange can be described as the fol- 



lowing 'three-point turn' sequence: 

ay-6 ^23^ y-6 ^12^ hy^ ^ bya (^^gj 

The sign change of 7a as it shifts from branch 2 to branch 
1 via branch 3 is given by the product of two signs: 

7a Saja , = S23S3I (50) 

For 75, the result is 

76^5676, Sb = Si2. (51) 

The arrow notation here implies that in computing the 
expectation value of a product of Majorana operators 
which includes 7^, 75 or both, the expectation value after 
the exchange will be multiplied by a sign equal to 5^,55 
or SaSb. Note that in our labeling scheme, the Majorana 
operators retain their original subscripts, even though 
their positions have been interchanged. We have now 
dropped the prime symbol used in the previous section 
to designate Majorana states in the final configuration. 

The values of the signs Sij can be obtained from the 
results of the previous section, using a permutation of 
the labels 1,2,3 where necessary. Each step in (49) is 
composed of two elementary parts. For example, in the 
first step, Majorana a is pulled back through the junc- 
tion from wire 2 to wire 1; then it is pushed through the 
junction onto wire 3. Thus we have 523 = S21S13. We 
may treat this process as a continuous rotation of the 
Majorana from the initial direction = (/>2 — tt to the 
final direction (/)3 — tt. The rotation must be done in the 
counterclockwise direction, so that aa is never equal to 

— TT. The sign S23 will be negative if a crosses the neg- 
ative X-axis in the process. Similarly, S12 will be negative 
if ai) crosses the negative x-axis when it is rotated in a 
counterclockwise sense from its initial direction (/)i — tt to 
(j)2 — TT. The sign 531 will be negative if crosses the 
negative x-axis as it is rotated from (/)3 — tt to — tt. If 
we choose the wire directions as suggested by the letter 
F, say with (01,02,^3) = (7r/6, 57r/6, 37r/2), we find that 

S23 = S12 = 1, S31 = -1, so Sa = -1, Sb = 1. 

Note that if we repeat the entire process twice, thereby 
restoring a and b to their original wires, both Majorana 
operators will be multiplied by -1. This follows from the 
fact that SaSb = —I. If we consider other orientations of 
the three wires, the sign Sa and may change, but their 
product will always be -1. 

In the second case of interest to us, the two Majoranas 
are on the edges of two different segments, which is de- 
scribed by the F-configuration ^Y^. [See panel (b) of 
Fig. 6.] Again, a three-step turn can exchange the two 
Majoranas: 

dyb ^23^ tyb ^12^ byt ^31^ byd ^^2^ 

In the first step, Majorana a is pushed through the junc- 
tion from wire 2 to wire 3; in the second step Majorana 
b is pushed from wire 1 into the junction and pulled onto 
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wire 2; and finally Major ana a is pulled from wire 3 to 
wire 1. The associated sign S23 will be negative if the 
angle passes through the negative x-axis while rotat- 
ing in a clockwise sense from the initial direction 02 to 
the direction (j)^ — tt. The sign S31 will be negative if 
aa passes through the negative x-axis when rotating in 
a clockwise sense from ^3 — tt to 0i. Finally, S12 will 
be negative if passes through the negative x-axis on 
rotating from (pi to ^2 • This rotation must also be taken 
in a counterclockwise sense, so that the avoids the 
direction ^2 — tt. 

As in the previous case, we may define signs Sa and 
55 characterizing the change in the Majorana operators 
after the transformation, with Sa = 523S31 and S5 = S12. 
Again, we find that the product SaSb is always negative. 
With the particular choice of orientations specified above, 
we have ^23 = s^i = 1 and S12 = —I 

As a direct application of the above results, we may 
consider what happens if we begin with a configuration 
like that shown in panel (b) of Fig. 6, and we are given 
the number parities P2 and Pi of two occupied segments 
in the initial state. If we perform two clockwise Majorana 
exchanges of the type described above, we will reach a fi- 
nal state in which the expectation values of P2 and Pi will 
each have changed sign. By contrast, if we perform only 
a single Majorana exchange. Pi and P2 will each have 
zero expectation value in the final state. The state will 
actually be a coherent superposition, with equal weights, 
of a state where the parities have changed and one where 
they are the same as in the initial state. 

As pointed out in Ref. 25, it is evident from our anal- 
ysis that the choice of chirality of the underlying p-wave 
proximity affects the exchange process and determines 
which Majorana obtains a minus sign upon exchange. In 
the N = 2 model, it is our choice for the sign of the 
spin-orbit coupling that governs this chirality. 



VI. SUMMARY 

The results of the previous analyses may be summa- 
rized as follows. For an arbitrary continuous manipu- 
lation of the parameters of the Hamiltonian, consistent 
with the requirement that the Majorana fermions always 
remain far apart, we may define a number qj for each 
Majorana, which equals ±1 depending on whether the 
parameters defining that Majorana have crossed an even 
or odd number of cuts where the Majorana wave function 
changes sign. (This must include any sign change that 
occurs when the Majorana is passed through a Y junc- 
tion). If the changes in the Hamiltonian are adiabatic in 
the sense that they are slow compared to the scale set by 
the gap to finite energy excitations, but fast compared 
to the exponentially small interactions between the sepa- 
rated Majoranas, then the expectation value of a product 
of any even number of Majorana operators 7^ taken af- 
ter the manipulation will be identical to the expectation 
value before the manipulation, except for a sign factor Q 



which is the product of the qj for the Majoranas involved. 

Suppose that at the end of a series of adiabatic ma- 
nipulations the system is divided into a set of disjoint 
TS wire segments. The quantities of immediate physical 
significance will be the electron number parities in the 
individual wire segments. The operator which measures 
the parity of a TS segment may be expressed as a product 
^7a76 of the Majorana operators at the two ends of the 
segment with an overall sign kab that depends on param- 
eters such as the orientations of the wire ends, but does 
not depend on the past history. Rules for computing the 
factors kab^ as well as the history-dependent Majorana 
factors qj have been given above, using a particular con- 
vention for the signs of the Majorana operators. If these 
factor are known, and if the expectation value (^7^7/0)0 
in the initial state is known for the two Majoranas which 
wind up, respectively, at the two ends a and b of the 
wire segment in question, we may calculate the expecta- 
tion value of the parity operator in the final state as the 
product of the factors kabQjQk and the initial expectation 
value {ijjjk)o' 

In the simplest case, where the ending Hamiltonian is 
the same as the starting Hamiltonian, so that the system 
contains the same geometric arrangement of TS clusters 
that it had initially, and if each Majorana is returned to 
the same position Rj that it had initially, the factor k 
will be the same in the initial and final states. Moreover, 
if the parity of each wire is known in the initial state, then 
the initial expectation values (i7j7/c)o may be deduced. 
The parity of a wire segment in the final state will then be 
equal to the parity in the initial state multiplied by the 
history-dependent product qjqk for the two Majoranas 
at the end of the segment. If the Majoranas at the ends 
of a wire segment originated on the same segment, but 
are reversed in position in the final state, the parity in 
the final state will still be determined by the parity in 
the initial state and the factor qjqk^ but there will be 
an additional factor of —1 due to the reordering of the 
operators. 

If the Majorana operators at the end of a wire seg- 
ment in the final state originated on two different wire 
segments in the initial state, the parity in the final state 
will not be determined by the parities of individual wires 
in the initial state. If the initial state had a definite 
parity on each wire segment, then the expectation value 
for the parity of the wire segment in question will be 
zero in the final state. However, the expectation value of 
the product of the parities of several such wire segments 
may still be non-zero (and equal to ±1), if the Majoranas 
at their ends, collectively, are a permutation of the Ma- 
joranas from an equal number of wire segments in the 
initial state. 

The analyses we have presented can be applied to the 
case of a three dimensional network of wires, as well as 
to a two-dimensional network. In all cases, the result 
of a set of adiabatic manipulations will be a topological 
invariant of the path in parameter space between the 
initial and final states. If the Majorana manipulations 



18 



are performed on a planar network of wires and junctions, 
with e perpendicular to the plane, with b parallel to e 
and with the phase of A the same for the entire system, 
and if the configuration of TS wires is the same in the 
initial and final states, then the result for any physical 
quantity will be the same as one would have obtained 
by analysis of the braiding of the positions Rj. These 
results would be the same as for the braiding of e/4-quasi- 
particles in non-Abelian v = b/2 fractional quantum Hall 
state^^ 
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